Search for: All records

Abstract We generalize a magnetogram-matching Biot–Savart law (BSl) from planar to spherical geometry. For a given coronal current densityJ, this law determines the magnetic field $\tilde{\mathit{B}}$whose radial component vanishes at the surface. The superposition of $\tilde{\mathit{B}}$with a potential field defined by a given surface radial field,B_r, provides the entire configuration whereB_rremains unchanged by the currents. Using this approach, we (1) upgrade our regularized BSls for constructing coronal magnetic flux ropes (MFRs) and (2) propose a new method for decomposing a measured photospheric magnetic field as $\mathit{B}={\mathit{B}}_{\mathrm{pot}}+{\mathit{B}}_T+{\mathit{B}}_{\tilde{S}}$, where the potential,B_pot, toroidal,B_T, and poloidal, ${\mathit{B}}_{\tilde{S}}$, fields are determined byB_r,J_r, and the surface divergence ofB–B_pot, respectively, all derived from magnetic data. OurB_Tis identical to the one in the alternative Gaussian decomposition by P. W. Schuck et al., whileB_potand ${\mathit{B}}_{\tilde{S}}$are different from their poloidal fields ${\mathit{B}}_{\mathrm{P}}^{<}$and ${\mathit{B}}_{\mathrm{P}}^{>}$, which arepotentialin the infinitesimal proximity to the upper and lower side of the surface, respectively. In contrast, our ${\mathit{B}}_{\tilde{S}}$has no such constraints and, asB_potandB_T, refers to thesameupper side of the surface. In spite of these differences, for a continuousJdistribution across the surface,B_potand ${\mathit{B}}_{\tilde{S}}$are linear combinations of ${\mathit{B}}_{\mathrm{P}}^{<}$and ${\mathit{B}}_{\mathrm{P}}^{>}$. We demonstrate that, similar to the Gaussian method, our decomposition allows one to identify the footprints and projected surface-location of MFRs in the solar corona, as well as the direction and connectivity of their currents. 

Abstract Global solar photospheric magnetic maps play a critical role in solar and heliospheric physics research. Routine magnetograph measurements of the field occur only along the Sun–Earth line, leaving the far side of the Sun unobserved. Surface flux transport (SFT) models attempt to mitigate this by modeling the surface evolution of the field. While such models have long been established in the community (with several releasing public full-Sun maps), none are open source. The Open-source Flux Transport (OFT) model seeks to fill this gap by providing an open and user-extensible SFT model that also builds on the knowledge of previous models with updated numerical and data acquisition/assimilation methods along with additional user-defined features. In this first of a series of papers on OFT, we introduce its computational core: the High-performance Flux Transport (HipFT) code (https://github.com/predsci/hipft). HipFT implements advection, diffusion, and data assimilation in a modular design that supports a variety of flow models and options. It can compute multiple realizations in a single run across model parameters to create ensembles of maps for uncertainty quantification and is high-performance through the use of multi-CPU and multi-GPU parallelism. HipFT is designed to enable users to write extensions easily, enhancing its flexibility and adaptability. We describe HipFT’s model features, validations of its numerical methods, performance of its parallel and GPU-accelerated code implementation, analysis/postprocessing options, and example use cases. 

The Sun’s corona is its tenuous outer atmosphere of hot plasma, which is difficult to observe. Most models of the corona extrapolate its magnetic field from that measured on the photosphere (the Sun’s optical surface) over a full 27-day solar rotational period, providing a time-stationary approximation. We present a model of the corona that evolves continuously in time, by assimilating photospheric magnetic field observations as they become available. This approach reproduces dynamical features that do not appear in time-stationary models. We used the model to predict coronal structure during the total solar eclipse of 8 April 2024 near the maximum of the solar activity cycle. There is better agreement between the model predictions and eclipse observations in coronal regions located above recently assimilated photospheric data. 

Abstract The Wang–Sheeley–Arge (WSA) model has been in use for decades and remains a popular, economical approach to modeling the solar coronal magnetic field and forecasting conditions in the inner heliosphere. Given its usefulness, it is unsurprising that a number of WSA implementations have been developed by various groups with different computational approaches. While the WSA magnetic field model has traditionally been calculated using a spherical harmonic expansion of the solar magnetic field, finite-difference potential field solutions can offer speed and/or accuracy advantages. However, the creation of new versions of WSA requires that we ensure the solutions from these new models are consistent with established versions and that we quantify for the user community to what degree and in what ways they differ. In this paper, we present side-by-side comparisons of WSA models produced using the traditional, spherical harmonic–based implementation developed by Wang, Sheeley, and Arge with WSA models produced using a recently open-sourced finite-difference code from the CORHEL modeling suite called POT3D. We present comparisons of the terminal solar wind speed and magnetic field at the outer boundaries of the models, weighing these against the variation of the WSA model in the presence of small perturbations in the computational procedure, parameters, and inputs. We also compare the footpoints of magnetic field lines traced from the outer boundaries and the locations of open field in the models. We find that the traced field-line footpoints show remarkable agreement, with the greatest differences near the magnetic neutral line and in the polar regions. 

Abstract Coronal Mass Ejections (CMEs) are immense eruptions of plasma and magnetic fields that are propelled outward from the Sun, sometimes with velocities greater than 2000 km/s. They are responsible for some of the most severe space weather at Earth, including geomagnetic storms and solar energetic particle (SEP) events. We have developed CORHEL-CME, an interactive tool that allows non-expert users to routinely model multiple CMEs in a realistic coronal and heliospheric environment. The tool features a web-based user interface that allows the user to select a time period of interest, and employs Regularized Biot-Savart Law (RBSL) flux ropes to create stable and unstable pre-eruptive configurations within a background global magnetic field. The properties of these configurations can first be explored in a zero-beta magnetohydrodynamic (MHD) model, followed by complete CME simulations in thermodynamic MHD, with propagation out to 1 AU. We describe design features of the interface and computations, including the innovations required to efficiently compute results on practical timescales with moderate computational resources. CORHEL-CME is now implemented at NASA's Community Coordinated Modeling Center (CCMC) using NASA Amazon Web Services (AWS). It will be available to the public in early 2024. 

Abstract One systematic limitation of solar coronal hole (CH) detection at extreme ultraviolet (EUV) wavelengths is the obscuration of dark regions of the corona by brighter structures along the line of sight. Another problem arises when using CHs to compute the Sun’s open magnetic flux, where surface measurements of the radial magnetic field, $B_r^{\odot }$, are situated slightly below the effective height of coronal EUV emission. In this paper, we explore these two limitations utilizing a thermodynamic magnetohydrodynamic (MHD) model of the corona for Carrington rotation (CR) 2101, where we generate CH detections from EUV 193 Å images of the corona forward-modeled from the MHD solution, and where the modeled open field is known. We demonstrate a method to combine EUV images into a full Sun map that helps alleviate CH obscuration called theminimum intensity diskmerge(MIDM). We also show the variation in measured open flux and CH area that is due to the effective height differences between EUV and $B_r^{\odot }$measurements. We then apply the MIDM method to SDO/AIA 193 Å observations from CR 2101, and conduct an analogous analysis. In this case, the MIDM method uses time-varying images, the effects of which are discussed. We show that overall, the MIDM method and an appreciation of the effective height mismatch provide a useful new way to extract a broader view of CHs, especially near the poles. In turn, they enable improved estimates of the open magnetic flux, and help facilitate comparisons between models and observations. 

Abstract We describe, test, and apply a technique to incorporate full-Sun, surface flux evolution into an MHD model of the global solar corona. Requiring only maps of the evolving surface flux, our method is similar to that of Lionello et al., but we introduce two ways to correct the electric field at the lower boundary to mitigate spurious currents. We verify the accuracy of our procedures by comparing to a reference simulation, driven with known flows and electric fields. We then present a thermodynamic MHD calculation lasting one solar rotation driven by maps from the magnetic flux evolution model of Schrijver & DeRosa. The dynamic, time-dependent nature of the model corona is illustrated by examining the evolution of the open flux boundaries and forward-modeled EUV emission, which evolve in response to surface flows and the emergence and cancellation flux. Although our main goal is to present the method, we briefly investigate the relevance of this evolution to properties of the slow solar wind, examining the mapping of dipped field lines to the topological signatures of the “S-Web” and comparing charge state ratios computed in the time-dependently driven run to a steady-state equivalent. Interestingly, we find that driving on its own does not significantly improve the charge state ratios, at least in this modest resolution run that injects minimal helicity. Still, many aspects of the time-dependently driven model cannot be captured with traditional steady-state methods, and such a technique may be particularly relevant for the next generation of solar wind and coronal mass ejection models. 

Abstract Differential emission measure (DEM) inversion methods use the brightness of a set of emission lines to infer the line-of-sight (LOS) distribution of the electron temperature (T_e) in the corona. DEM inversions have been traditionally performed with collisionally excited lines at wavelengths in the extreme ultraviolet and X-ray. However, such emission is difficult to observe beyond the inner corona (1.5R_⊙), particularly in coronal holes. Given the importance of theT_edistribution in the corona for exploring the viability of different heating processes, we introduce an analog of the DEM specifically for radiatively excited coronal emission lines, such as those observed during total solar eclipses (TSEs) and with coronagraphs. This radiative-DEM (R-DEM) inversion utilizes visible and infrared emission lines that are excited by photospheric radiation out to at least 3R_⊙. Specifically, we use the Fex(637 nm), Fexi(789 nm), and Fexiv(530 nm) coronal emission lines observed during the 2019 July 2 TSE near solar minimum. We find that, despite a largeT_espread in the inner corona, the distribution converges to an almost isothermal yet bimodal distribution beyond 1.4R_⊙, withT_eranging from 1.1 to 1.4 in coronal holes and from 1.4 to 1.65 MK in quiescent streamers. Application of the R-DEM inversion to the Predictive Science Inc. magnetohydrodynamic simulation for the 2019 eclipse validates the R-DEM method and yields a similar LOST_edistribution to the eclipse data. 

Abstract We present in this Letter the first global comparison between traditional line-tied steady-state magnetohydrodynamic models and a new, fully time-dependent thermodynamic magnetohydrodynamic simulation of the global corona. To approximate surface magnetic field distributions and magnitudes around solar minimum, we use the Lockheed Evolving Surface-Flux Assimilation Model to obtain input maps that incorporate flux emergence and surface flows over a full solar rotation, including differential rotation and meridional flows. Each time step evolves the previous state of the plasma with a new magnetic field input boundary condition, mimicking photospheric driving on the Sun. We find that this method produces a qualitatively different corona compared to steady-state models. The magnetic energy levels are higher in the time-dependent model, and coronal holes evolve more along the following edge than they do in steady-state models. Coronal changes, as illustrated with forward-modeled emission maps, evolve on longer timescales with time-dependent driving. We discuss implications for active and quiet Sun scenarios, solar wind formation, and widely used steady-state assumptions like potential field source surface calculations. 

Heavy ion signatures of coronal mass ejections (CMEs) indicate that rapid and strong heating takes place during the eruption and early stages of propagation. However, the nature of the heating that produces the highly ionized charge states often observed in situ is not fully constrained. An MHD simulation of the Bastille Day CME serves as a test bed to examine the origin and conditions of the formation of heavy ions evolving within the CME in connection with those observed during its passage at L1. In particular, we investigate the bimodal nature of the Fe charge state distribution, which is a quintessential heavy ion signature of CME substructure, as well as the source of the highly ionized plasma. We find that the main heating experienced by the tracked plasma structures linked to the ion signatures examined is due to field-aligned thermal conduction via shocked plasma at the CME front. Moreover, the bimodal Fe distributions can be generated through significant heating and rapid cooling of prominence material. However, although significant heating was achieved, the highest ionization stages of Fe ions observed in situ were not reproduced. In addition, the carbon and oxygen charge state distributions were not well replicated owing to anomalous heavy ion dropouts observed throughout the ejecta. Overall, the results indicate that additional ionization is needed to match observation. An important driver of ionization could come from suprathermal electrons, such as those produced via Fermi acceleration during reconnection, suggesting that the process is critical to the development and extended heating of extreme CME eruptions, like the Bastille Day CME.